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Electromagnetically induced transparency (EIT) is an important tool for controlling light propa- 
gation and nonlinear wave mixing in atomic gases with potential applications ranging from quantum 
computing to table top tests of general relativity. Here we consider EIT in an atomic Bose-Einstein 
Condensate (BEC) trapped in a double well potential. A weak probe laser propagates through one 
of the wells and interacts with atoms in a three-level A configuration. The well through which the 
probe propagates is dressed by a strong control laser with Rabi frequency S7 M , as in standard EIT 
systems. Tunneling between the wells at the frequency g provides a coherent coupling between 
identical electronic states in the two wells, which leads to the formation of inter-well dressed states. 
The macroscopic interwell coherence of the BEC wave function results in the formation of two ultra- 
narrow absorption resonances for the probe field that are inside of the ordinary EIT transparency 
window. We show that these new resonances can be interpreted in terms of the inter-well dressed 
states and the formation of a novel type of dark state involving the control laser and the inter-well 
tunneling. To either side of these ultra-narrow resonances there is normal dispersion with very 
large slope controlled by g. For realistic values of g, the large slope of this dispersion yields group 
velocities for the probe field that are two orders of magnitude slower than standard EIT systems. 
We discuss prospects for observing these ultra-narrow resonances and the corresponding regions of 
high dispersion experimentally. 



I. INTRODUCTION 

Electromagnetically induced transparency (EIT) [l[ is 
a quantum interference effect that occurs in coherently 
prepared three-level A atomic systems. The great utility 
of EIT comes from the fact that an ordinarily opaque 
medium can be made transparent to a probe laser while 
at the same time having large controllable dispersion and 
large third order nonlinear susceptibilities 0, @]- EIT 
uses a strong control beam to dress an electronic excited 
state with a third auxilliary level. A weak probe field, 
which normally has only a single excitation path from 
the ground state to the excited state in the absence of 
the control beam, now has two excitation pathways to 
the excited state, corresponding to the two dressed states 
formed with the auxilliary state. The resulting destruc- 
tive quantum interference between excitation pathways 
leads to vanishing absorption at the bare atomic reso- 
nance. Along with the vanishing of the probe absorption, 
the real part of the linear susceptibility, 5i[x ]> exhibits 
normal dispersion with a very large slope leading to ex- 
tremely slow group velocities for the probe field 0, 
Slow light propagation through EIT systems has been 
observed experimentally in a variety of media, including 
hot atomic gases [6[ and atomic Bose-Einstein conden- 
sates (BEC's) [7|, and is now well understood. 

Slow light propagation in EIT can be thought of in 
terms of the quasiparticles known as "dark states polari- 
tons" [8, 9], which are a superposition of the probe pulse 



and the atomic polarization of the ground states. Be- 
yond the novelty of simply controlling the speed of light, 
EIT has found numerous potential applications. Dark 
state polaritons provide a method for fully reversible stor- 
age of light pulses in an atomic medium by adiabatically 
switching on and off the control laser TC|. Light storage 
has important applications for quantum information pro- 
cessing since quantum information can now be transmit- 
ted by flying qubits (photons) between stationary qubits 
(atomic ensembles) in a quantum network. Besides engi- 
neering the linear susceptibility, EIT results in construc- 
tive quantum interference for the nonlinear susceptibil- 
ity, x^ 3 \ m the middle of the transparency window where 
the absorption vanishes and the dispersion is large. Such 
large nonlinearities in lossless media lead to an efficient 
scheme for four-wave mixing and frequency conversion 
in atomic vapors [l], HH, 03 • Additional work has shown 
that these large nonlinearities can be used to achieve non- 
linear mixing between pulses involving a few photons [l3| , 
which could be used to create an all-optical controlled- 
NOT gate 14], the essential element of a quantum com- 
puter. Furthermore, one of the most striking applications 
of EIT has been the realization that the propagation of 
ultraslow light in moving atomic media is mathemati cally 
the same as light propagation in curved spacetime [15j . 
Leonhardt and Piwnicki [l6| showed that in this case a 
vortex, such as in an atomic Bose-Einstein condensate, 
will behave like a black hole for the light. This opens up 
the possibility of table-top tests of general relativity. 
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In the current paper, we describe a modification of the 
standard 3-level EIT configuration that utilizes coherent 
tunneling of a BEC in a double well potential and leads 
to qualitative changes in the linear susceptibility of the 
probe laser, which, as a result, provides additional con- 
trol over the dispersion. More specifically, we consider 
the optical properties of an atomic BEC [17] with three 
electronic states in a A-configuration that is trapped in 
a double well potential [H, [H, [2(J HH ■ One well is pre- 
pared as a standard EIT system: the electronic excited 
state is coupled to one of the two stable ground states 
via a strong control laser, while a weak probe couples 
the other ground state to the same excited state. Both 
lasers are confined to a single well, leaving the second well 
unperturbed by them. However, the barrier between the 
wells represents a weak link through which atoms can 
tunnel. The global phase coherence of the condensate 
wave function can lead to phase coherent tunneling of the 
condensate wave function between the wells. This tun- 
neling is the origin of Josephson oscillations of the pop- 
ulation difference between the wells, which have recently 
been observed in a double well condensate [HI, [HI [2(| ■ 

This double well BEC 'Josephson junction' signifi- 
cantly modifies the probe EIT spectrum since the tunnel- 
ing transforms the 3-level A system of a lone well into a 6- 
level system spatially distributed between the two wells. 
The additional eigenstates of the 6-level system mani- 
fest themselves in the form of new absorption resonances 
in the probe susceptibility. In particular, we show here 
that the tunneling induces two qualitatively new ultra- 
narrow absorption resonances situated in the middle the 
EIT transparency window with widths and positions de- 
termined by the tunneling frequencies between the wells. 
These new resonances are clearly visible when the tunnel- 
ing frequencies are much less than the control laser Rabi 
frequency such that the transparency window is larger 
than the separation between these resonances. At the 
same time, the change in the probe index of refraction 
to either side of these new resonances is more dramatic 
still than found in a standard EIT system with a control 
laser of the same intensity, providing additional control 
over the group velocity and shape of a light pulse trav- 
eling through the medium. We predict that for realistic 
tunneling rates (~ lKHz), the dispersion to either side 
of the new resonance can yield group velocities up to two 
orders of magnitude slower than would be possible in EIT 
without the tunneling. Although one can in principle re- 
duce the group velocity in EIT by simply turning down 
the intensity of the control laser, this approach narrows 
the transparency window and is fundamentally limited by 
ground state decoherence. By contrast, the ultranarrow 
tunneling resonances presented here are independent of 
the ground state decoherence and can be used to create 
ultra-slow group velocities independently of the control 
laser strength. 

The remainder of this paper is organized as follows. 
In section II, we will describe our model for the two-well 
A BEC dressed by both the control beam and tunnel 



coupling, which has an analytic steady-state solution to 
its master equation. In section III, we will derive the 
system's linear susceptibility, x > from which we can 
extract the absorption coefficient and dispersion. In sec- 
tion [iVl we will consider the prospects for experimental 
observation of the ultra-narrow features we describe, and 
their consequences for slow light. 

II. THE MODEL 

The present contribution concerns a gas of N weakly 
interacting atoms of a Bose-Einstein condensate trapped 
in two neighboring wells of a double well potential. There 
have been various experimental realizations of double 
well potentials for BEC's involving some combination of 
magnetic and/or optical dipole potentials. The first of 
these used a focused blue-detuned far-off resonant laser 
in the center a harmonic magnetic trap [221 ] . Later at- 
tempts created double well potentials via two parallel 
laser beams that generated adjacent optical dipole traps 
within the same condensate [23. l24j . In these cases, tun- 
neling between wells was negligible. More recently, a 
double well potential with coherent quantum mechani- 
cal tunneling of the condensate wave function between 
wells was demonstrated [HI, [ill [5(| ■ These represented 
the first realizations of a single Josephson junction in an 
atomic BEC and serve as a guide for our EIT model, 
since coherent coupling of the wells is the essential new 
element. In these experiments [19j, the double well was 
created by superimposing a one dimensional optical lat- 
tice on top of the harmonic optical dipole trap leading to 
a potential in the x-direction, 

V e (x) = Knuj e x 2 + V t cos 2 (jj^J (1) 

where di is the lattice constant and £ is the electronic 
state of the atoms. This is because, in general, any mag- 
netic or optical potentials used to trap the atoms will 
depend on their electronic state and therefore atoms in 
different states will experience slightly different trapping 
potentials. We assume that in the z-direction, the har- 
monic trapping potential is much weaker than in the x 
or y directions, leading to an elongated cigar shaped po- 
tentials for the two wells with the long axis along the 
z-direction. 

We consider three internal electronic states of the 
atoms in a A configuration, denoted by eigenkets \a), \b), 
and |c) where \a) is an electronically excited state, while 
\b) and |c) are hyperfine ground states of the atoms. The 
direct transition between the two lower levels is assumed 
to be dipole forbidden while the transition between the 
highest level and each of the lower levels are allowed opti- 
cal dipole transitions. Here we use a ' to denote the same 
internal states but in the left well so that for example, 
| a) is the electronic excited state in the right well while 
| a') represents the same internal state of the atom but 
now that atom is located in the left well (see Fig. [TJ. 
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FIG. 1: Schematic description of our system: the atoms in 
the right well are dressed by a strong control beam (indicated 
here by its frequency, tJ M ) near resonance with the \c) — » |o) 
transition. We are interested in the behavior of a weak probe 
beam (here, uj p ) propagating in the right well, near resonance 
with the b) — > \a) transition. Atoms in electronic state \£) 
are coupled via tunneling through the inter- well barrier to the 
corresponding states \£') in the left well. The lasers propagate 
along the axis perpendicular to the page. 



In analogy to the standard EIT configuration, we as- 
sume one of the two wells (in this case the right well) 
is dressed with a strong control beam with electric field 
amplitude £^ and frequency lo^ that is close to resonance 
with the energy difference between levels \a) and |c). 
Here we are concerned with the propagation through the 
right well of a weak probe field, £ p , with frequency uj p 
near resonance with the |6) — > |o) transition. 

The restriction that the lasers interact with only a sin- 
gle well should be achievable provided the spacing be- 
tween the wells is sufficiently larger than the diffraction 
limit. The diffraction limit is essentially given by the 
wavelength of the probe and control lasers, which we de- 
note simply as A. Based on Eq. [TJ the well spacing must 
satisfy d > A, which can be achieved with current tech- 
nology. For example, in the experiment of Ref. [l8| the 
spacing between the wells is 4.4/xto, which is significantly 
larger than a typical optical wavelength. Furthermore, 
the group in Ref. [20| were able to optically resolve a 
single well to successfully image tunneling effects. The 
probe and control lasers are assumed to propagate along 
the z-axis to maximize the optical thickness of the sam- 
ple. 

As we are working in the zero-temperature limit, we 
adopt the Hartree approximation and assume that all of 
the particles are co-existent in a single fully-condensed 
state. We model the wells as weakly coupled har- 
monic potentials [2l[ with ground state wave functions, 

u[ L ^ R \r), localized in the left (L) or right (R) wells, 
which also depend on the electronic state since they rep- 
resent the localized ground state near the minima of the 
state dependent potential, Vg. We assume that the over- 



all condensate wave function, ^(r, t) can be expressed in 
terms of these basis functions ug(r), 

9(r,t)=VN E Mt)4 R) (*)\Z)+ E Mt)tf\r)\t) 



£—a.b,c 



' —a'b'c' 



As a matter of notation, we introduce the vectors 
of probability amplitudes for the right and left wells, 
respectively, ^ = (ipa,ipb,ipc) T and likewise W = 

We work in a rotating frame defined by 



where 



R = 







(3) 



In this basis, the Gross-Pitaevskii equations for the six 
probability amplitudes are: 
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(4b) 



(4c) 



(4d) 



(4e) 
(4f) 



We have assumed that these amplitudes are normalized 
to 1: 

E \M*)\ 2 + E \M*)\ 2 = i 



Here we have incorporated the ground state ener- 
gies of atoms in the wells, J d 3 r[w^ fe ''(r)]*[— ?i 2 V 2 /2m + 

Vi(v)]v!(^ (r), into the definition of the atomic energy lev- 
els, uj£. The couplings between levels are moderated by 
their complex Rabi frequencies defined as Ml ac e~ t ^ ac — 
£)j,D ac for the control field and Ml ab e~ l<t ' ab — £ p D ab 
for the probe field. Here f2y is taken to be real and 
Djj = e (i \ x • e \j) is the dipole moment matrix element 
in the direction of the laser polarization, e. 
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In principle, each of the atomic levels is subject 
to a different coupling constant for the tunneling be- 
tween wells. We denote these couplings by hgi/ 2 = 
- J d 3 r[u[ L \r)]*[-h 2 V 2 /2m + Vtir^u^M. For the 
sake of completeness, we note that Ref. [25| shows that 
hgi is equal to the Josephson coupling energy, Ej, that 
appears in the Hamiltonian for the bosonic Josephson 
junction (2(1 ]. The two-body interactions are denoted by 
a rank 3 tensor, U^fe, defined by: 










" 




u, = 










(5) 















where the index i runs likewise over a, 6, and c. The el- 
ements U i:j = (AnhaijN/m) f d 3 r\u\ k) (r)\ 2 \uf ] (r)\ 2 rep- 
resent the interaction strength (in units s _1 ) between 
states \i) and \j) in the same well in terms of the s-wave 
scattering length between the two states, ay. 

The absorptive and dispersive properties of the 
medium with respect to the probe are given by the linear 
susceptibility, which we can derive from the coher- 
ence between states |o) and \b) . To proceed further we 
must introduce the density matrix p, defined as the outer 
product of the probability amplitudes, 



(6) 



For the sake of clarity we note that ('J, ty') = 
(ip a ,ip b ,ip c ,^ a ',ip bl ,i>c')- By the product rule, we arrive 
at the rate of change of the density matrix: 

dt dt { ' ' 1 ' ' dt 2 1 ,Ph 

To incorporate decay, we have introduced the decay ma- 
trix r whereby each element of the density matrix decays 
at the rate p tj cx -JijPij where 7y = (7, +7j)/ 2 + l[^ P) ■ 
Here, 74 for i = a, 6, c, a' ' ,b' ,c' are the decay rates for 
the populations and 7^- is decoherence due to pure de- 
phasing for i j. Extending our notation, we write the 
density matrix in the rotating frame defined in equation 
Eq. [3] as p. Extracting the equation of motion for p abl 
we find: 



1 dt 



= {\-i lab + ¥{\J a -XJ b )^)p ab 

+ — e (Paa - Phh) - — e ^""Pch 



9b - 



9a - 
2 Pa'b, 



(8) 



where we have defined the probe's detuning from the 



j b) transition, A f 



U>b 



Before proceeding further, we note that equation © 
depends on five additional, mutually dependent linked 
differential equations, each of which is likewise coupled 
to other terms in the density matrix. But we are only 
interested in the linear susceptibility for the probe field, 



so we may solve the resulting coupled system of equations 
to first order in the strength of the probe field £ p , which 
we have assumed to be weak. We assume that initially all 
of the atoms are in the two states \b) and \b') . As a result, 
up to order £%, we have p aa = p a , a , = p aa , = p a > a = 0. 
Additionally, |c) only develops population at order £p£ 2 
in perturbation theory and therefore to first order in the 
probe laser, p cc = p dd = p ac = p a , c = p ad = p a , c , = 
Pcd = Pdc = 0. The only terms in the density matrix 
that are nonzero to zeroth order in the probe are p bbl 
Pdb', Pb' b and pbb' while to first order in the probe p ab , 
Pa'b, Pa'd, Pad, Pcb, Pdb, Pcd, and p dV are nonzero. 

The control laser is assumed to be of arbitrary strength 
so that we must solve the equations to all orders in Q ac . 
In addition to this, we solve to all orders in the tunneling 
rates gi. This is because the critical element of EIT is 
the presence of coherence between the two ground states, 
p cb . In the case that the tunnel coupled states \b) and 
\b') as well as |c) and |c') are nearly degenerate, they will 
form superposition states between the two wells that will 
in turn affect p cb . It is important to point out that our 
choice to include g b and g c to all orders is not at odds 
with our choice to only keep terms to linear order in the 
probe despite the fact that g b and g c are themselves small. 
The assumption of a weak probe means that p aa <C p bb 
at all times. To second order in perturbation theory one 
can easily show that p aa ~ (£! a fc/7a) 2 - Tunneling re- 
sults in finite populations for both wells and in the case 
of degenerate states (including mean field interactions), 
there is equal population in both wells p bb — p b > b > = 1/2. 
Consequently, as long as there is finite population in \b), 
the weak probe condition remains Q ab <C 7a and is only 
weakly effected by gi. Appendix B gives the full solution 
for Paa to second order in fl ab . 

In order to keep the inter-well couplings to all or- 
ders, we move to a partially dressed state basis, in which 
the {\b) , \b')} and {\c) , |c')} subspaces of our effective 
Hamiltonian are diagonalized. To simplify matters, we 
take x dp ^ — 76 = 75/ = j c — -f d — 0, which is a reason- 
able approximation because the decay rates for atoms in 
a BEC are given by the lifetime of the condensate, which 
is much longer than all other time scales in this prob- 
lem. We keep the decay from the excited electronic state, 
7a = 7a' 7^ 0, which is due to spontaneous emission. 

The effective Hamiltonian for the {\b) ,\b')} subspace 
can be written as a sum of its diagonal and traceless 
parts: 



H 



1,1/ 



h [ U! b 

h 
2 



^' ] \j b ^' 



-9b 



~9b 
, T ~,t, 



(9) 



where 8w — ^U&'I' — \J b ^' is the energy difference 
between the corresponding states in each well. In the 
case that p bb = p b ' b >, 5 bb i = 0. If the wells are initially 
prepared with equal population in both of them, then 
S b d = initially and will remain zero since the eigen- 
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states of Hbb' have equal probability to be \b) and \b') in 
this case. 

The matrix diagonalizing Hbb' will be a member of 
SO (2), which can be written in terms of a rotation angle 
in the {\b) ,\b')} subspace: 



where 



n _ f cos 9b sin0 o x 
6 — sin 9b cos 9b j 



COS b = 9? y /2 
sin0b = (i±^Z£rV /2 

„eff 



(10) 



9? - VtiV 



(11a) 

(lib) 
(11c) 



Under this transformation, we find dressed states \B) 
and \B') whose probability amplitudes can be written in 
terms of the bare states and the angle of rotation 9 b : 

4> B = cos 9 b ipb + sin 6 Vv (12a) 
■4> B > = - sin 9 b ipb + cos 9 h i> b > (12b) 

Identical reasoning applies for the {|c) ,|c')} subspace 
leading to the dressed states {|C) , \C) }. Note that these 
dressed states are in a coherent super-position of spatially 
delocalized states. 

Combining these transformations, we arrive 

at the full transformation from the bare basis 

{\a) , 1 6) , |c) , \a') , 1 6') , \c') } to the dressed basis 
{\a) ,\B) ,\C) ,\a') ,\B') ,\C')}, 



D = 











cos 9b 



— sin 9b 
\0 






cos 9. 



— sin 9 C 



\ 

sin(9 6 

sinfl c 

1 
cos6» b 

cos 9 C J 



(13) 



The Gross-Pitaevskii equations rewritten in the dressed 
basis are given in appendix A. 

The transformation for the density matrix from the 
original basis, p, to the dressed basis is pd = DpD^ . In 
terms of the dressed states, the coherence p a b is given by 



p ab = COS 9 b PaB ~ Sm0 b p aB '. 



(14) 



Our assumption that all atoms are initially in some com- 
bination of { 1 6) ,!&')}, implies that in the dressed basis 
the terms pbb, Pb'B', and pb'b = (pBB')* are in general 
nonzero to zeroth order in the probe. Beginning with 
these and keeping only terms to first order in fl a b, we 



arrive at two decoupled systems of four equations each: 



■ d PaB f A , gf 1 



dt 



A p + -Ubb + Uab - hub ) PaB 

Qac —i(h„ r I n ~ ■ n ~ \ 9a ~ 

- ~7T e (cos 9 c p C B - sm9 c pc'B) - y/VB 



(cos9 b p~BB - sin9 b p B 'B) (15a) 



„eff „eff 2 1 



n 



(15b) 



.dp, 



C'B 



dt 



.dp a > 



^e**« cos9 c p aB 

off off 1 1 

A p - A, + ^- + g -f + -U cb - -Ubb ) Pcb 



+ ^e^sin^ aB 



dt 



ff 1, 



Q 

A p + -| ^U b b + U ab - hab ) Pa'B 

9a - 
~^PaB\ 



(15c) 



(15d) 



and likewise 

.9 paB 



. dpCB' 

1 dt 



= (a p 


gf 


2 


- — e 




Qab 
+ ~ e 




= Up 


-A,- 



1 



-U bb + U ab - ij a b PaB 



9a ~ 



off off 1 i 



(16a) 



(16b) 



.dpa'B' ( . gf l TT , TT . \ - 

I gt = I Ap — - — 1/ 66 + £/„(, - «7 afc I pa'B' 



5a - 

y/VB'- 



(16c) 



(16d) 



where the control laser detuning is = Lu a — uj c — lu^ . At 
this point we assume that the zeroth order populations 
in the dressed states are nonzero and controlled by a 
tunable parameter, ip, such that p BB — cos 2 (6 b — ip) and 

P^B'B' = sm2 (^6 — f)- At the same time, we assume 
that coherences between the dressed states are initially 
zero, Pg) B = Pbb' = 0- This is a reasonable assumption 
since if the atoms are specifically prepared at some time 
in the past in the dressed states or simply allowed to 
equilibrate to the eigenstates of the double well, then any 
coherences would be destroyed before the experiment by 
even a small amount of decoherence. The effect of initial 
coherences between dressed states on the transient probe 
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absorption spectrum in a three level system has been 
considered before [27j and shown to give rise to temporal 
oscillations in the absorption coefficient similar to optical 
nutation. 

The structures of the solutions to these are identical. 
In both cases, we write the systems as linear equations 
of the form dX/dt — — M • X(t) + A, and note that such 
equations have steady-state solutions of limt_>oo X(t) = 
M _1 • A. The necessary terms p a B and p a B' are then 
the corresponding elements of the resulting vectors. The 
general form of the analytic solution is quite complicated 
and is given in appendix B. 



Likewise, the complex linear susceptibility is given by 
— V/^eoSp), which determines both the absorption 
coefficient, a(uj p ) — kpSslx^ 1 '] and the index of refraction, 



n(ujp) 



3?[x (1) ]) 1/2 . The spatial term [u a (r)]*u b (r) 



reflects the density profile of the atoms and only deter- 
mines the optical thickness of the condensate. 

The simplest case to consider is when the dressed state 
mixing angles are 9 b = 9 C = 7r/4, corresponding to sym- 
metric and antisymmetric superpositions between the 
two wells. This occurs when 8w = &cc' — or equiva- 
lently when p bb — pvv ■ In this case the solution simplifies 
considerably and the polarization is given by 



III. OPTICAL PROPERTIES OF THE RIGHT 
WELL IN THE DEGENERATE ENERGY CASE 

The polarization for the probe is related to p a b by 

V = 2N[u a {v)]*u b {v)D ab ~ Pab . (17) 



N[u a {v)]*u b {v)Dl b £ p 



z. 



(18) 



where 



Z± = 



((2A M - 2A P ± g b f - g 2 c ) (2A P T 9b - 2i lab )(l =p sin2<^) 



((2A„ - 2A p ± g b f - gl) ((2A p T 9b - 2i lab ) 2 



9l) 



(2A M - 2A P ± <? 6 )(2A P T 9b — 2ij ab )Q 



Note that here we have redefined the probe and control 
laser detunings to include the mean field energy shifts, 
A p + U a b Ap and A M + U ab — > Ap. We emphasize the 
fact that, in the limit that g a , g b , <? c , f — > 0, we recover 



the standard EIT coherence, which for A^ 
form 



has the 



Pab 



2 (Ap(Ap - i lab ) - (n ac /2) 2 ) 



Further analysis requires us to estimate the values for 
the important variables in the problem. Atomic sponta- 
neous emission rates for atoms commonly used in BEC 
experiments (Na, Rb, Li) are typically on the order of 
lOMhz. Typical tunneling times were on the order of 
10ms in Ref. [l8| while more recent experiments achieved 
tunnel couplings between the wells as high 7900Hz (26|. 
We therefore assume that lKHz is a reasonable estimate 
for the coupling strength between two wells of a BEC. 
Therefore, unless stated otherwise, we will use the values 
gi = 10 3 s _1 and j a = 10 7 s _1 , which yields gi = 10~ 4 7 a . 
Examples of the real and imaginary parts of the suscep- 
tibility are shown in Figs. 2 and 3 for A M = while 
Figs. 0] and [5] display the spectrum's dependence on the 
tunneling parameters. (Note that in these and all sub- 
sequent figures the susceptibility is plotted in units of 
N[u a (r)]*u b (r)D^ b /2e h'y ab .) We see that the presence 
of the second well manifests itself as two ultranarrow res- 
onances located inside of the EIT transparency window. 
When S bb i — 0, the new resonances are symmetrically lo- 



cated about Ap = at the locations ±g b /2. In general 
the location of these resonances is A p = ±<?b /2 and for 
f! ac , "f ab 3> g b -,g c their shape is approximately Lorentzian 
with a full width at half max of 



9c 

&>ac 



lab 



(19) 



Note that we have also obtained this same result even 
when we have included decoherence due to pure dephas- 
ing between the ground states \b) <-> \b'), \c) «-> |c'), and 
| b) <-> |c) [28]. This implies that these resonances will 
be clearly separated even in the presence of finite ground 
state decoherence, and, in particular, 7^. Similar results 
have been obtained by Lukin et al. [291 ] and Mahmoudi 
et al. who studied EIT in 4-level system where an 
additional ground state was coupled via an RF or optical 
transition to the same ground state that is coupled to the 
control laser. They found an ultra-narrow resonance in 
the EIT spectrum with a line width of the same form as 
Eq. [TO] 

These new resonances can be understood in terms of 
the interaction of the dressed states of the \b) and \ b') sub- 
system with the eigenstates of the |c') <-> |c) <-> \a) sub- 
system. To first order in the probe, \a') can be neglected 
altogether. First let us consider how the tunnel coupling 
would effect the probe absorption for the \b) — * \a) tran- 
sition without any control laser. The dressed states \B) 
and \B') are both populated ground states that couple 
directly to the excited state by the probe laser. The ener- 



FIG. 2: The full EIT spectrum of the system near the \a) -\b) 
resonance. The 3[x ] is plotted as a dotted line, and 5R[x' 1 ' ) ] 
as a solid line. Note the two additional features, symmetri- 
cally located around zero detuning, at ±<?(,/2, with equal am- 
plitude for if — 0. In this plot we have taken gt — g c = fab/W 
(~ bOOKHz) to emphasize the modifications to the standard 
EIT spectrum. Given more physically realistic parameters, 
the features would be considerably narrower and closer to- 
gether. Cf. Figs. [3] |4] and [5] Note that here and in subse- 
quent figures Q. ac = j a = 2-y ab . 




FIG. 3: Close-up of the positive detuning narrow absorption 
peak corresponding to the presence of the second well. As 
in the previous figure, S^ 1 '] is plotted as a dotted line, 
and Jftfx' 1 '] as a solid line. Note that to either side of the 
absorption peak, the real part of the linear susceptibility is 
rapidly changing, whereas the absorption goes to zero still 
more quickly. In this plot, gb — g c = {^ab) x 10 -4 , which 
we estimate as a reasonable upper bound for the coupling 
between wells (see text). Here ip = 0. 

gies of \B) and \B') are Hujb.b' — ^>b^9h /2. Therefore 
even in the absence of the control laser, the |b) — > \a) ab- 
sorption line would be split into two new lines located 
at uj a — cob and to a — u)b>- This is essentially an Autlcr- 
Townes doublet induced by the tunneling. Figure[6]shows 
the two ultranarrow resonances as a function of <p, which 
controls the relative population in the dressed states such 
that for <p = it/ 4, pg B = 1 while for ip = 3-7T/4, Pg) B , = 1. 

The excited state \a) is in fact coupled to \c) via the 
control laser while |c) is coupled to \d) via the tunneling. 



y ab 0.0001001 



FIG. 4: Ssfx' 1 ^! which is proportional to the probe absorp- 
tion coefficient, plotted near the positive detuning resonance. 
Here <p — 0, gt = (S'fab) x I0 -4 is fixed, and g c is varied, 
to demonstrate how g c moderates the width of the tunneling 
induced resonances. See Fig. 0for a cross-section of this plot 
at g c = lKHz. 
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-0.0002 -0.0001 0.0001 0.0002 

Yab 

FIG. 5: S[x (1) ] but now gb is varied, keeping g c = (27 ab ) x 
I0 -4 fixed along with tp = 0. gb moderates the distance be- 
tween the peaks. As gb goes to zero, the peaks merge to create 
a single narrow peak at the origin. 



This system is a three-level system that is isomorphic 
to a A atom. Again assuming 6 cc i = and the control 
laser is on resonance, A M = 0, then we have the following 



0.0002 
1.0 




FIG. 6: Q[x W ] for 9b = 9a = {2~1ab) x 1CT 4 as a function 
of ((3. One can see that the amplitude of each resonance is 
proportional to the initial population in \B) and \B'). 



Hamiltonian for the {|a), |c), |c')} subsystem, 



H = hw n l- 



o n ac o 

-,g c 



(20) 



The eigenstates of this Hamiltonian are 



1 



= —= (sin0|a) + \c) + cos 9\c'}) 
V2 

= -= (sin0|a) - |c) + cos 0\c')) 
v2 



\ a o) 



where tanfl 



cos u\a) — sin B e 



(21) 

(22) 
(23) 



-n„ 



s/gc. The energies of the states \a±) 
are E± — hui a ± fty'f^ + while |ao) has energy 
Eq = fku a . As one can see, |ao) is the same type of dark 
state that appears in STIRAP and coherent population 
trapping. In this case, this tunneling induced dark state 
is a superposition of \a) and |c') but not |c) . Since it is 
decoupled from the control laser, there will not be any 
destructive quantum interference in the probe absorption 
for transitions to \ao). Transitions from the {\B), \B')} 
manifold to |ao) will then exhibit absorption resonances 
at oj a — u>b,b'> which correspond to the new ultranarrow 
resonances. This independence of |ao) from |c) also ex- 
plain why the line width, r„, does not depend on either 
led or 7b c . This is in stark contrast to what would hap- 
pen if g c = 0, which would correspond to a conventional 
EIT system but with two ground states, {\B), \B')}. In 
this case, destructive interference created by the control 
would lead to nulls in the absorption at uj a — ojb.b 1 - 

A general understanding of locations of the absorp- 
tion resonances can be obtained from Fig. [7] which 
shows a schematic diagram of the energy levels of the 
dressed ground state manifold {\B), \B')}, which are cou- 
pled to all three states of the excited state manifold 
{|ct+), |<Z-)j | fl o)} v i a the probe. All in all there are six 



| Bo). E «=f"»a 



h i 

| a.), E_=ha>~Jg] + 



\B), E B =hm b +hg b l2 
\b), E b =ha t 

\B), E B =hm b -Ks b l1 



FIG. 7: Energy level diagram that indicates transitions in- 
duced by the probe laser between the ground state manifold 
{\B}, \B')} and the excited state manifold {|a+}, \a~), |ao)}- 
Transitions to the dark state \ao) are indicated by dashed 
lines. The energy of the bare state \b) is also shown for refer- 



transitions that should appear as resonances in the ab- 
sorption spectrum. The transitions to the two 'bright' 
states \a±) correspond to the main absorption peaks lo- 
cated at A p ss ±f2 ac /2 for £l ac ^> g Cl gb- Notice that 
each of these resonances actually consist of a pair of res- 
onances separated by a distance gb but because 7 a h 3> gb 
these pairs cannot be individually resolved. Since to first 
order in the probe field \a) is never populated, the tun- 
neling rate between \a) and \a'), g a has only a negligible 
effect on x*' 1 ''- 

As we can see from Eq. [18] there is a transparency 
window of width g^ in between the two ultranarrow res- 
onances with Q : [x^ 1 ' ) ] = at A p = 0. To cither side 
of these resonances, the absorption is negligible. In the 



vicinity of these resonances the dispersion, d^t[x^']/ 9u 
is extremely large and, to either side of the peak, there 
is a region of width w g c /2 in which the dispersion is 10 
times greater than in standard EIT while within a region 
of width f» <7 c /8, the dispersion is amplified by a factor of 
100. (Note that 'Standard EIT' refers to the case where 
.9a = gb = 9c — but with all other parameters, including 
the control laser, being the same.) The absorption, mean- 
while, drops to below 1% within an order of magnitude 
of the feature's width — 2(<7 c /f2 ac ) 2 7 a {, — from the center 
of the resonance, and is negligible within the regions of 
interest. Using our approximation of g c ~ lKHz, we find 
there is a region, to either side of the peaks, of width 
O(lKHz) in which the absorption is negligible (< .001%) 
and the dispersion is ~ 10 times greater than in standard 
EIT with a control laser of the same strength; similarly 
there is a region of width O(100Hz) in which the dis- 
persion is « 100 times greater than standard EIT, again 
with negligible absorption (likewise < .001%). This im- 
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0.000105 



0.00011 
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FIG. 8: The group velocity in a window of 100Hz just to the 
side of one of the ultranarrow resonances, in units of the stan- 
dard minimum EIT group velocity for a given atomic species 
and control beam strength. Here gb = g c = (27 a b) X 10 -4 and 
(/3 = 0. As can be seen, the group velocity is suppressed by 
a factor on the order of 10 2 , relative to standard EIT. The 
group velocity of standard EIT, weit, is calculated in the ab- 
sence of any tunnel coupling g a — gt — g c = 0. Absorption is 
negligible in this frequency window (see Fig. [3J) 



plies that the group velocity 0, 0] , 

c 



n+(uj p /2n)(dn[x^}/du; p ) 



could be made significantly smaller than in previous ex- 
periments (see Fig. [8]). Note that the slope of the disper- 
sion is rapidly changing within each of these regions, and 
so any pulse transmitted through the well would undergo 
considerable reshaping. 



IV. DISCUSSION AND CONCLUSIONS 

The regions in which dispersion is especially high in 
this system are on the order of g c or smaller, begging the 
question of their experimental accessibility using readily 
available equipment. We propose two solutions. The first 
follows a recent paper by Pfau's group, at Stuttgart [3ll |. 
In their experiment, only a a + polarized probe couples 
the initial state of the atoms to the same excited state 
as the control laser. On the other hand, a a~ polarized 
probe couples off-resonantly to other excited states that 
are unaffected by the control laser. By mixing a small 
amount of o~ polarized light into their otherwise er + po- 
larized probe beam, they simultaneously measured the 
absorption and dispersion of the EIT system by examin- 
ing the interference pattern between the a + and er - po- 
larized components of the probe. They report observing 
features as narrow as 4KHz in the Sft^ 1 )]. Performing a 
similar experiment with a system prepared as described 
in this paper would test our predictions for the modified 
absorption and dispersion arising from the presence of 
the second well. 



To directly observe the low group velocity that follows 
from our predictions would require lasers with line widths 
small relative to the frequency window in which the dis- 
persion is large. Though such lasers are not readily avail- 
able commercially, several groups have reported perform- 
ing spectroscopic experiments within the range of inter- 
est. As early as 1999, for instance, Young et al. at NIST 
were able to achieve sub-hertz width lasers, albeit with 
nontrivial active stabilization [32| . More recently, using 
a a 657-nm diode laser with a femtosecond comb, a col- 
laboration at NIST and LANL performed KHz-resolution 
spectroscopy on cold neutral calcium (33|. And so, the 
next generation of lasers could take full advantage of the 
high-dispersion regime exhibited by our double well sys- 
tem. 

Finally we would like to comment on the approxima- 
tions made in our model of a double well BEC that could 
affect the experimental feasibility of our proposal. The 
major assumption made is that the spatial profile of the 
condensate is fixed. This means that the ground state 
wave functions u[ L ^ R \r) are unchanging and also that 
there are no excitations of the condensate. Similar 'two- 
mode' models for a double well BEC have been success- 
fully used to obtain quantitative agreement with experi- 
ment [n m, m . 

In the limit of small nonlinear interactions between 
atoms, u[ L ' R '(r) are the single particle wave functions, 
Gaussians in the case of our harmonic potential. The 
shape of these wave functions are not affected by the 
population in the wells. However, in the limit of large 
nonlinear interactions, uf 1 (r) are best approximated 
by Thomas-Fermi wave functions and therefore will only 
remain the same if the number of atoms in each well docs 
not change with time. This is implicit in our choice of 
initial conditions and a weak probe field that does not 
significantly excite \a). Furthermore, the dynamics of \a) 
are dictated by the laser coupling, which is faster than 
the time needed for atoms in |o) to equilibrate in the 
potential. Atoms excited from \b) to \a) will not have 
time to equilibrate in the new potential and therefore 
^^^'(r) = ui L (r) implying that the prefactor in Eq. 
[T71 is maximal: [w a (r)]*Uf,(r) — + |wf,(r)| 2 . 

The second issue is that of excitations of the con- 
densate. Transverse excitations can easily be avoided 
by using tight confinement in the x and y directions 
(uj x ,ujy > 10 3 s -1 ). However, since we want to maximize 
the optical thickness of the wells along the z direction, we 
must have very weak confinement along this axis, imply- 
ing that excitations along this direction could readily be 
excited. These long wavelength excitations should have 
a negligible impact for two reasons. First, the measured 
signal is proportional to the integral of the condensate 
density in the z-direction. Even if there are local spatial 
variations of the density in the z-direction, those varia- 
tions would be averaged out in the integration. Secondly, 
the trapping frequency along the z-axis will be at least 
one to two orders of magnitude smaller than in the x and 
y directions (oj z ~ 10s _1 — 100s _1 ). The momentum dis- 
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tribution of such elongated quasi-lD condensates have 
been measured experimentally using Bragg spectroscopy 
[35l ] and showed to have a momentum distribution in the 
axial direction ranging from 50Hz to 500-Hz for tem- 
peratures ranging from T = 0.257c to T = 0.9Tc and 
trapping potentials lo x = lo v = 4700s -1 and uj z = 30s" 1 
(T<7 is the critical temperature for Bose-Einstein conden- 
sation). Since experiments have already been performed 
with interwell tunnel couplings as high as 7900Hz [2^ |. 
this would imply that at low temperatures the separation 
between the ultranarrow resonances would be more than 
100 times greater than the inhomogeneous broadening in 
the z-direction. Additionally we can point to the work in 
Ref. that measured the affect of thermal excitations 
in a double well BEC Josephson junction. This work 
showed that the loss of inter-well phase coherence due to 
thermal excitations was negligible when the tunnel cou- 
pling energy was much larger than the thermal energy, 
fcgT. Therefore we conclude that for T -C To, inhomo- 
geneous broadening due excitations will be sufficiently 
small as to not affect our results. 

In conclusion, we have studied the effect of interwell 
tunneling on the EIT dispersion and absorption of a 
probe laser interacting with a double well Bose-Einstein 
condensate. As a future direction we plan to extend these 
results to an optical lattice where the ultranarrow reso- 
nances should be determined by the band structure of 
the lattice. Additionally, the effect of tunneling on the 
nonlinear susceptibility will be explored. 



.dtpa' 

.dips' 

l ~dT 



.dtpc 

l -d7 



(uj a - W/t - w p + #' f U a #) - y^a (24d) 



• • ^(* tU fc* + *"U**')J ipB' 

(24e) 



^e*-' sin Ma 



off 1 

LJ C - UJ p + ^- + - (*t Uc * + & 1 U c *0 ) i>c> 



xTV'^TT \T>'" 



^e**" sin Ma- 



(24f) 



The rotation angle (9& and g e h ^ are defined in the text. 
The transformation to the dressed state basis {|C), \C')} 
from {|c) , |c') } is defined in the same manner as Eq. 
(10) with the rotation angles and dressed state energies 
explicitly given by 



cos 6 r 



i-s~s/g? y /2 

off \ 1/2 



9f= V^+al 



(25a) 

(25b) 
(25c) 



APPENDIX A 



Here we first present the Gross-Pitaevskii equations in 
terms of the dressed states of the subspaces {\b) ,\b')} 
and {\c) , \d) } 



and S cc i = 



(* + U« 



* - ^^Uc*') /2. 



- ^e~ i0 » 6 (cos 9 b ipB - sin Mb 

- ^e" 4 ^ (cos Mo - Bin Mo 
ff 1,= 



.dip 1 



dt 



.dtpc 



(24a) 



fib 



(* t U 6 * + *' t U 6 ^') ) %pB 



(24b) 



^£ e **« COS Ma 



(24c) 



VI. APPENDIX B 



The general solution for the coherence p a b is given by 
the steady state solution of Eqs. (JTSJ) and (TIB)) combined 
with Eq. (|14p. The general form of the steady state p a b 
is then: 
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Pab = 



^ ( C ° s2( ^"^ )CQs20fc ((2A, - 2A P - 9 f f - (gff) (2A P + gf - 2 llab ) 



sin 2 (9 b V ) sin 2 9 b ({2 ^ _ 2 ^ + gf) , _ {2Ap _ gf _ 2iiab) j (26) 



C+ 

while the population in \a) to second order in the probe field is 
iQ 2 ab f f cos 2 (6b — if) cos 2 9b 



P. 



a a 



7a V V C- 

sin 2 (6*b - ip) sin 2 b 



((2A„ - 2A p - .gf ) 2 - (gf ) 2 ) (2A p + .gf - 2z 7ab ) 



((2A M - 2A p + gf) 2 - (gf) 2 ) (2A p - gf ~ 2 ?7ab ) 
(2A M - 2A p - gf) 2 - (gf) 2 ) (2A p + gf + 2i lab 



C+ 

cos 2 (9 b -f) cos 2 6> b f /oA oA cffN 2 
C* 

+ Sin2( ^ Sln2 ^ ((2A M - 2A p + .gf) 2 ~ (9f) 2 ) (2^-^ + 2***))) (27) 
Here we have defined 

C± = ((2A„ - 2A p ± gf ) 2 - (gf) 2 ) ((2A p T 5 f - 2 ?7ab ) 2 - 5 2 ) 

+ (2A p T fff - 2i 7 ab) (2A M - 2A p ± 5h - g c cos(29 c )) Q 2 ac (28) 

We have simplified equations ([2^]) through by making the substitutions A p — ^ + [7 ab — * A p and A p — + ?7 ab — > 
A,- 
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